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Abstract 

Using the fact that the c-quark inside a charmonium moves with a small velocity v in the charmonium 
rest-frame, one can employ an expansion in v to study decays of charmonia and results at the leading order 
for XcO, 2 77 exist in the literature. We study corrections at the next-to-leading order in the framework of 
nonrelativistic QCD(NRQCD) factorization. The study presented here is different than previous approaches 
where Xc0,2 is taken as a bound-state of a cc pair and a nonrelativistic wave-function is used for the pair. We 
find that the corrections are consist not only of relativistic corrections, but also of corrections from Fock state 
components of Xc0,2 in which the cc pair is in a color-octet state. For Xc2 there is also a contribution from a Fock 
state component in which the pair is in a F-wave state. We determine the factorization formula for decay widths 
in the form of NRQCD matrix elements representing nonperturbative effects related to Xc0,2, and calculate the 
perturbative coefflcients at tree-level. Because the NRQCD matrix elements are unknown, a detailed prediction 
for the decay Xc0,2 — > 77 can not be made, but the effect of these corrections can be determined at certain level. 
Estimations show that the effect is significant and can not be neglected. 



PACS numbers: 13.25.Gv, 14.40. Gx, 12.38.Bx 
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Table 1: List of major experimental results 









CLEO(1994) 

CLEO(2001) 

E835(2000) 

L3(1999) 

OPAL(1998) 


3.76 ± 1.07 
1.61 ± 0.75 


1.08 ± 0.28 
0.53 ± 0.16 
0.270 ± 0.042 
1.02 ± 0.25 
1.76 ± 0.37 



1. Introduction 



Because of their simple final states, decays of P-wave charmonium into two photons can be good channels to 
determine the value of ag and for study of properties of charmonium. In recent years various experiment groups 
have pubhshed their results on Tjy{xco) and Tj^{xc2) 0, |[ ||, Some major and newest results are hsted in 

Table 1, where errors from different sources are combined. From these data it is easy to find that although there is 
a tendency of consistency in latest results of CLEO and E835, values of different groups have a large discrepancy 
from each other. On the theoretical side, because charm quark can be taken as heavy quark and it moves with a 
small velocity v inside a charmonium in its rest-frame, one may describe the c- or c quark inside charmonia with 
a nonrelativistic wave-functions by taking charmonia as a bound-state of a cc pair. In most previous calculations 
for decays of charmonium such a nonrelativistic wave- functions is employed for Xc0,2 ~^ 77[^ @i i[ ^ 0- By 
expanding the small velocity v, the decay width at the leading order of v can be expressed with the derivatives of 
the wave functions at origin, where the leading order is at u^. The corrections from the next-to-leading order of 
V are also studied with the nonrelativistic wave-functions iQ, H, ^, in which the corrections are only relativistic 
corrections. This also implies that at the next-to-leading order the charmonia Xc0,2 are still taken as a bound-state 
of a cc pair and the pair has the same quantum number of XcO.2- 

From point of view of a quantum field theory like QCD, a hadron state is a bound state of quarks and gluons 
and has many components composed with different number of quarks and of gluons. For charmonia it can be 
expanded in v and for XcO.2 it looks like: 

|XcO,2) = l(cc)i) + 0{v)\{cc)sg) + higher orders, (1) 

where the index 1 or 8 denotes color singlet or color octet respectively. From this expansion the state |(cc)i) 
is a dominant component, while other components have a suppressed probability to be observed in Xc0.2i the 
suppression parameter is v. Hence, in a systematic study of higher-order corrections one should also include 
contributions from these suppressed components, i.e., a systematic expansion in v is needed. Recently, a novel 
approach based on NRQCD for productions and decays of quarkonia is proposed jl2|. In this approach a decay 
width can be systematically expanded in the velocity v. For example, the decay width of XcO can be written as: 

^{XcO ^ U) =^Cn{Xco\On\Xco), (2) 
n 

where On's are operators defined in NRQCD and Cn can be calculated with perturbative theory. There is a rule 
of power counting of for the operators 0„ |ll| . We will follow this approach to study the next-to- leading order 
corrections for the decay XcQ.2 77- We will determine the form of relevant operators in NRQCD and calculate the 
corresponding coefficients. It should be noted that the role of the suppressed component |(cc)83) has been studied 
before in inclusive processes. It has been shown |l3| that only by including the contribution from \(cc)gg) the 
inclusive decay of Xcj can be consistently predicted with QCD. It has been also shown that such components with 
a color octet cc pair of ■0 are very important for explaining the t/j surplus at the Tevatronp^. In this work we will 
show the effect of such components in exclusive processes. The effect of \{cc)8g) has been studied in Xcj — > Trvrp^], 
in which a light-cone wave-function for the component was introduced and the effect was not parameterized with 
matrix elements of NRQCD. 
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Figure 1: Diagrams of contribution for decay width from |(cc)i) state 

At the next-to-leading order of v, the corrections, as we will see, consist of relativistic corrections and also 
corrections from the state \{cc)sg)- Beside these, for Xc2 there is also a correction from the state |(cc)i) where 
the cc pair is in a state with the orbit angular momentum I = 3. These corrections take a factorized form as 
products of pertubative coefficients and matrix elements, in which the matrix elements are defined in NRQCD and 
they parameterize nonperturbative effects. We will determine the form of these matrix element and calculate these 
coefficients at leading order of as- Our work is organized as the following: In Sect. 2 we calculate these coefficients 
and give our main results. In Sect. 3 we discuss the effect of these corrections and Sect. 4 is our summary. 



2. The correction at the next-to leading order 



We consider XcO.2 to two photon decay in its rest frame: 

XcoMP) ^ lih) + 7(fc2) (3) 
where P, /ci, k2 denote momentum of Xc0,2 and two photons. We decompose the decay width as 

r^Ti + rg, (4) 

where Fi g denotes the contribution from a color-singlet- and color-octet cc pair, respectively. The color-singlet 
contribution can be written as: 

Ti = 16 QMJ dV J / d^xe^^^^^+^^y^ j d^ye'^^^('^-yHm{^)cAv)\Xc0.2) 



■pi- rric pr- h- "^c ^ pi - ^t- mc y^, 

where c{x) stands for the Dirac field of c-quark and i, j, i' and j' stand for color- and spin indices. This contribution 
can be represented by the diagrams in Fig.l. The matrix element with Dirac fields of c-quark represents the 
nonperturbative effect related to the initial hadron. The quantity H is calculated with perturbative theory. An 
expansion in v can be performed by expanding the Dirac fields with fields of NRQCD. This can be done by using 
Foldy-Wouthuysen transformation||l^ : 
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(8) 

where ^-nd x(2^) ^'I'e NRQCD fields and ipix') annihilates a c(c)-quark respectively, rric is the pole mass of 
c-quark. Using this transformation the matrix elements are then of the NRQCD fields x, and their derivatives. 
Operators D and di are of order and Dq and do are of v"^ order in NRQCD according the power counting rulep2|. 
The space-time dependence of matrix element are controlled with different scales in the different directions. In 
the time direction it is rricv'^, while in the spacial direction it is rricV. The expansion in v is then completed by 
doing Taylor expansion of fields x "0 around origin. Using the power counting rule in [l^ we collect all local 
matrix elements up to order of w^, which are allowed by rotation-, charge conjugation- and party symmetries for 
Xc0,2, and calculate the corresponding contributions. Doing the calculation in this way is similar to the matching 
method by taking a free cc pair, as proposed in the difference is that we only collect the contributions from 
operators with correct quantum numbers in the problem considered here. The calculation is rather tedious but 
straightforward. We obtain: 



ri(xco^77) = ^^^(Xco|Oi^A./fPo)|Xco)-^^^^^(Xco|WrPo)|Xco) + 0(«'^) 



■S \X.cO\yEM(^ ' ^"^^-^^ 

rric 

(9) 

Oem{'Po) = i^^(^D.^)x|0)(0|xt(^D.,T)V (10) 

GemCPo) = l[l^HY^nY^-'r)x\0){0\xHY^-'r)^ + h.c.], (11) 

where (Ji{i — 1,2,3) is the Pauli matrix. With the above results we find the correction to the color-singlet 
contribution is only the relativistic correction. Similarly, the decay width of Xc2 is: 

ri(Xc2 > ll) = ^^^^^(Xc2|C£;Af(^^2)|Xc2) " ^"^^^^^^ (Xc2 I^BM (^^2) |Xc2) 

- 'Tj^sty (Xc2\g'EMeP2)\Xc2) + 0(.^) (12) 

OEMeP2) = Vt(_!B(v^-))x|0)(0|xt(^B(V))V' (13) 
WfP2) = i[^t(_!e)2(_!S(V.))^|0)(o|xt(^S(V))^-t-/i.c.] (14) 

where the notation (ij) means 

Tiv) {T'^ +TJ')/2 - (5yT'='=/3, i.e., the tensor is symmetric and trace-less. Unlike 
the case with Xco we find not only a relativistic correction but also a contribution from a mixing with a _F-wave cc 
pair. This is allowed by symmetries, the cc pair is in a spin triplet state, a spin triplet state with a orbit angular 
momentum I = 3 can have a total angular momentum J = 2, 3, 4. The state with J = 2 is just the Xc2 state. 

Now we consider color octet contributions from |(cc)8.g). At the order we consider the contributions can be 
represented by diagrams in Fig. 2. The contributions before the expansion in v can be written as: 

AFs, = Qe^e^l rfr j d*y J d4^^e^'=fe-)+»(fci+fe^)-(0|c-,(x)c, (2/)|xcO,2) 

d'y'd'x'J ^^e^'='^'+^^-'(xcO,2|Q'(0)5.Gnx')c,.(y')|0) x M^,,^,,, (16) 



3 



with 



ATsb = Qe'e'J dV J d^y J d4^^e^'=(J/-)+^(fei+'==)-(0|c-i(:r)c, (2/)|xc0.2) 

X I d^y'd^x' j ^^e-^'='^'-^''^'(xcO,2|c7Ky')5«G''(a;Oc,'(0)|0) x M^,-,„,v (17) 
AFs, = Qe'e^l dV J d^y J d4^^e''=(J/-)+^('=i+'=^)-(0|c-,(a;)c, (y)|xcO,2) 



(19) 



1 1 

''''' ~ ft-jt^m, + ie^ ft-ft,-m, + ie^ \j 

I 1 1 , 

M'm'f = [7" , „ , , 7^ + 7^ 



(20) 



^ - It2- rric + ie ' ' ' - ^i- rUc + ie'^ 



(21) 



Similarly we need to perform the expansion in w for matrix element like {xco\ci{x)gsG^{y)cj(0)\0). Again with 
the power counting rule, one can find that the leading term for the matrix element is related to the matrix element 
of the operator ip^gg'E ■ crx- This matrix clement is at order of and will lead to a contribution to Fg at order of 
u^. To calculate the contribution we take the gauge: G° = 0, in this gauge the electric chromofield is: 

E=|g. (22) 

Performing the expansion for matrix elements in Eq.(19), results can be obtained in a straightforward way. We 
obtain the contribution at v"^: 

r8(Xc0^77) = =^^^^^^{Xco\:FEMePo)\Xco)+0{v'^) (23) 

:FemCPo) = i[V'V-s,Ex|0)(0|xV-DV' + /i.c.] (24) 

r8(Xc2— 77) = 5 {Xc2\TemCP2)\Xc2) + 0{v') (25) 

:Fem{^P2) = ^[VV(*5«E^)x|0)(0|xV(^Bj)V + /i-c.]. (26) 

Eq.(9), Eq.(12), Eq.(23) and Eq.(25) are our main results. From these results we can see that the corrections 

consist not only of relativistic correction which is represented by the matrix elements {xco\Gem{^ P())\Xc()) and 
{Xc2\Gem{^P2)\Xc2), but also of contributions from the color octet component. For Xc2 there is also a contribution 
from the component in which the cc pair is in color single state with the orbit angular momentum 1 = 3. 
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Figure 2: Diagrams contribute to decay width from \{cc)gg) state, the diagrams for conjugated contributions are 
not shown. 

3. Numerical impact of the correction 

To make theoretical predictions for decay widthes one needs to know all matrix elements in the last section. 
Unfortunately, there is no information available for the matrix elements appearing in the correction. There are 
estimations for the matrix element at the leading order of v, but for consistency we need an estimation for them at 
the accuracy of 0{v'^). All this prevents us from detailed predictions. However, numerical impact of the correction 
can be studied by noting that there are several relations among these operators. Because of spin symmetry of 
NRQCD we can have: 

(Xc2|Obm('P2)|Xc2) = {Xco\OEM{'Po)\Xco){l + Oiv')), 
(Xc2|W('P2)|Xc2) = (Xc0|W('P0)|Xc0)(l + O(«')), 

{Xc2\TEMeP2)\Xc2) = {xco\^emCPo)\Xco){1 + 0{v^)). (27) 

As discussed before, the first relation is not useful here unless we know the correction at 0(u^). Beside these 
relations one can also use equation of motion to obtain another type of relations. Generally there are subtleties for 
using equation of motion for operators. However, it is shown that one can still use equation of motion for operators 



which are sandwiched between physical states 18 . For operators relevant to J/^* one has used this to derive this 
type of relations [p7|. 

The equation of motion of NRQCD reads: 

(zA + ^)V' = 0, (iDt-^)x^O. (28) 

Zm Zm 

We can use these equations to replace the operator with Df. For doing this we note that the matrix element 



5 



appearing in Eq.(lO) can be written as: 

where the corrections terms are due to interchange of the gauge covariant derivative D, and they are at order of 
0{v'^) and can be neglected here. Using the equation of motion we obtain: 

(Olx^D • aD^V; + (B\)^B ■ <j^ijj\xco) = -2m,z[(0|x^D • aA^ + (Ax)^D • a^\xco)] 

= -2m,£;^^„(0|x^D-aV'|xco)-2m,(0|x^gE-aV'|xco), (30) 

where E^^^ is the binding energy. Doing the same for Xc2 we obtain the relations: 

(Xcol^('^'o)IXco) = rncE^^„{xco\OePo)\Xco)+ni,{xco\H''Po)\Xco), 

{Xcjg{^P2)\Xc2) = m,E^^,{XcJOCP2)\Xc2)+ma{Xc,\H''P2)\Xc2), (31) 

and the binding energy is defined as: 

Exao = Wxco - 2mc, £^^02 = Wxc2 - 2mc. (32) 

It should be noted that the difference E^^^ ~ P'xoo is at order higher than w^, hence the above relations are in 
consistent with those in Eq.(27). To study the numerical impact we define the following parameters: 

_ (Xc0|^e/'0)|Xc0) _ {Xc2\H^P2)\Xc2) (l + 0(„2)) 



m,(Xc0|O(3P0)|Xc0) me(Xc2|0(3P2)|Xc2) 

(Xc2|g'eP2)|Xc2) , . 

^ m,2(x,2|0(3P2)|Xc2)' ^ ^ 

These parameters are order of . With these parameters and with the relations in Eq.(31) the decay widths can 
be written: 

r(XcO— 77) = ^^^(xco|0(^Po)|Xco)-(l + 0.1787^^- 2.333:^ -3a8l+0(««), 



r(Xc2 77) = , % (Xc2|g(^P2)|Xc2) • 1 - 5.333^^^ - 1.024^ - 2.2738a8 - 1.6905aF 

+0{v^), (34) 

where we have added the one-loop correction to the contribution at leading order of If we take rric = 1.5GeV 
and as(TOc) ~ 0.3, we obtain by using experimental values of hadron masses: 

r(Xco^77) « ^^^(xcolOfPo)IXco)- {1 + 0.017 -0.645 - 0.3} + 0(^;6), 

r(xc2^77) « ^^^^4 (Xc2|0('P2)|Xc2) ■ {1 - 0.510 - 0.380 - 0.227 - 0.169} + 0{v^), (35) 

where we have taken = ap — Q.l. In the above equations, the numbers in brackets correspond to the terms in 
brackets in Eq.(34), respectively. From these numbers, we see that the effect of binding energies is very significant, 
also the effects from color octet contributions and from mixing with the F-wave state are not small, although their 
effects can not be determined precisely. However, our results show that each correction at the next-to-leading order 
is significant. Unless among those corrections some cancellation happens, they can not be neglected. 



4. Summary and conclusion 
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We have presented a systematic study of high order corrections in two photon decays of XcO,2- Using NRQCD 
velocity power counting rule, implicit expansion of decay width in relative velocity v is achieved. We calculated not 
only first order relativistic corrections but also corrections from suppressed component \{cc)sg) in which cc pair is 
in a color octet state. The effect of color octet has been studied extensively in inclusive processes. In this work we 
have studied the effect in exclusive processes. Unlike most previous theoretical calculations in which nonrelativistic 
wave function is employed, our calculation is based on a field theoretic framework of NRQCD. The results are 
expressed by products of perturbatively calcidable short distance coefficients and nonperturbative matrix elements 
which concern with initial bound states. In the case of Xc2, a correction from orbital angular momentum Z = 3 is 
also obtained in our calculation. In previous calculations only relativistic corrections can be estimated. 

Unfortunately, the matrix elements at the next-to-leading order are not known at all, although some relations 
among them can be obtained by symmetries and by using equation of motion. This prevents us from numerical 
predictions for the decays. However, numerical impact of these corrections can be estimated at certain level and 
with this estimation the effect of the corrections is significant and can be not neglected. It would be interesting 
to study these matrix elements by nonperturbative methods, especially by lattice QCD. With a nonperturbative 
method one may determine how important the color octet component is. 
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